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We discuss heating and decoherence in traps for ions and neutral particles close to metal- 
lic surfaces. We focus on simple trap geometries and compute noise spectra of thermally 
excited electromagnetic fields. If the trap is located in the near field of the substrate, the 
field fluctuations are largely increased compared to the level of the blackbody field, leading 
to much shorter coherence and life times of the trapped atoms. The corresponding time 
constants are computed for ion traps and magnetic traps. Analytical estimates for the size 
dependence of the noise spectrum are given. We finally discuss prospects for the coherent 
transport of matter waves in integrated surface waveguides. 
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1 Introduction 

In the field of particle cooling and trapping, a strong trend towards miniaturisation and integra- 
tion has emerged in the last few years. Small particle traps might form the building blocks of 
future quantum computers. Integrated atom optical circuits might distribute coherently matter 
waves for interferometric or nanolithographic applications. Major issues in this field are heat- 
ing and decoherence that have to be controlled in order to maintain the coherence properties as 
long as possible. Heating is an intriguing concern because of the dramatic temperature differ- 
ence between the trapped particles and the macroscopic objects that form a typical miniature 
trap. In fact, recent experiments with small ion traps have revealed that the life time of the ion's 
vibrational ground state gets shorter in smaller traps, making it difficult to down scale the trap 
geometry to micrometre size [|J g]. In neutral particle traps, sizes in the micrometre range have 
already been achieved experimentally [^|, [|, |(J|, but to our knowledge, life times are difficult 
to measure and are rather limited by background pressure and inelastic few-body collisions. 

In this contribution, we consider a simplified geometry for particle traps and compute heat- 
ing and decoherence rates. The work is divided in two parts. We start with the heating of the 
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vibrational motion of an ion in a tightly confining potential. Using perturbation theory, the 
heating rate is linked to the cross correlation spectrum of the electric field at the trap center. 
This noise spectrum is evaluated asymptotically, taking into account that the trap distance is 
small compared to the photon wave length associated with the ion's vibration frequency. As a 
by-product, we also get the noise spectrum for the magnetic field. This spectrum determines 
the life time of a neutral particle in a magnetic trap, since the fluctuating field induces spin 
transitions to a non-trapped state, kicking the particle out of the trap. In the second part, we 
focus on the quasi-free motion of a particle in a wave guide (linear or planar). The particle 
scatters from thermal field fluctuations and thus loses its spatial coherence. A transport theory 
is formulated and analytically solved in the limit of a broad-band fluctuation spectrum. 



2 Heating of a trapped ion 
2.1 Heating rate 

Let us focus on a single degree of freedom of the ion's motion and assume a parabolic confining 
potential. The Hamiltonian is then simply given by 

fftrap = m (tfb + |) (1) 

with the ion's vibration frequency (typically in the MHzjange). If the ion is perturbed by a 
time-dependent force, this is described by the potential 



V(t) = x • F(r, t) = -an • F(r, t) (6 f + b) . (2) 

We assume that the force varies on a spatial scale much larger than the size a = (h/ (MQ)) 1 / 2 
of the trap ground state and evaluate it at the trap center r. x = in is the ion's displacement 
from the center. 

Using standard second-order perturbation theory, one may easily derive a master equation 
for the reduced density matrix of the ion, that describes its dynamics when the fluctuations 
of the force field F(r,i) are traced over [Q, |§J ||. This master equation allows to derive 
the equation of motion for the population of the lowest trap levels, as well as for the average 
creation operator and the average excitation number (in the absence of cooling processes) 

Poo = -7-P00 + 7+P11 (3) 
(n) = -( 7+ -7_)(n) +7- (4) 

(b) = -iO(6)-i( 7+ - 7 _)(6) (5) 

The rates 7 ± in these formulas are related to the cross spectral density of the force fluctations 
as follows 

7± = 7(r; ±fi) 

,2 



^}_>n^(r;±0) (6) 



1 -j 
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with 

/+oo 
dr<F i (r,t + r)F i (r ! t))e i ^ (7) 
-oo 

According to (^J), the heating rate from the trap ground state is thus given by To— >i(r) = j_. 
We are thus left with the calculation of the spectral density (^|) for the perturbing force. 

2.2 Electric field noise spectrum 

The trapped ion being charged, it is perturbed by fluctuating electric fields, F = qE. The 
Planck blackbody spectrum then gives the following spectral density 



„2 * 3 Sij 



where T is temperature. One must keep in mind, however, that this gives the thermal spectrum 
only in free space, far from the sources. But the ion is trapped at a distance z from the trap 
electrodes that is typically much smaller than the photon wave length A = 2nc/Vl associated 
with the vibration frequency. It is thus located in the near field of the electrodes, and the Planck 
formula does not cover this case. This was recognised already in the early days of ion trapping 
when 'hot' ion clouds (with temperatures much larger than room temperature) were cooled 
down by thermalisation with the surrounding electrodes, the coupling being provided by the 
absorption of the electric fields radiated by the moving ions in the lossy metallic environment 
[pj. When laser cooling took over to reach temperatures in the /LtK range, voltage fluctuations 
due to electric losses became a source of heating. Modeling the ion trap as a lumped circuit 
with resistance R(lu), the standard Nyquist formula for Johnson noise gives an electric field 
spectrum 

q 2 k B TR(u) 

S F (r;uj) « (9) 

z z 

where the high-temperature limit ksT ^> U\oj\ has been assumed. For a given trap geometry, it 
is not easy to determine the resistance R(uj) that enters this formula. Assuming that the electric 
currents propagate only in the skin layer of the electrodes, the NIST group estimated that the 
spectrum (^|) actually gives a heating rate too small to account for the experimental observations. 
In addition, there are indications that the scaling law To-a (r) oc Sf (r) oc 1/ z 2 is not followed 
experimentally (although it is difficult to exclude other influences when down-scaling the trap 
geometry) [[|. 

We now outline a microscopic effective model |[ ||] that may in principle allow to 
compute the electric noise spectrum for an arbitrary trap geometry. We describe the electric 
properties of the surrounding metal by its frequency-dependent complex dielectric function 
e(r;u>). The microscopic source of Johnson noise are fluctuating polarisation fields P(r, t) 
that are thermally excited in the metal. According to the fluctuation dissipation theorem, the 
spectral density of this fluctuating polarisation is related to the imaginary part of the dielectric 
function @ [IJ, 0]: 
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Fig. Fig. 1. Sketch of the microscopic effective model used to compute the electric field noise spectrum. 



Maxwell's equations now determine the field radiated by this polarisation. It may be written as 
an integral over the Green tensor 



E i {v;u>)= /dr'^G^ry^P^a;) (11) 



where G 1 ^ depends on the geometry of the trap electrodes.^ 

To proceed with the calculation, we now fix the geometry and consider a trap located at a 
distance z above an infinite flat metallic surface. In this geometry, the Green function in ( |TT| ) is 
explicitly known in spatial Fourier space, and it is possible to perform an asymptotic expansion 
in the near field limit z <C A. As a result, we obtain the following interpolation formula ||[ ^ 

s e ( z ^)~ -r-r, — hu i Q n TN s (s ij + 5 ij -^-] . (12) 

As expected, the spectrum only depends on the distance from the surface. The material proper- 
ties enter through the specific resistance g and the skin depth 8(uj) — c^£qq/lj. The geometry 
enters through the power law Se{z) oc 1/z 3 in the 'extreme near field' z <C <5. This behaviour 
may be understood from the principle of detailed balance: the heating rate of the ion is equal to 
the relaxation rate of an oscillating electric dipole, and it is well known that close to a metallic 
surface, this rate is dominated by nonradiative transfer and increases as 1/z 3 Jl5| jig ]. 

We may also extract from ( |l2| ) an 'effective resistance': if the distance is large compared 
to the skin depth, the near field spectrum indeed follows a 1/z 2 power law, as in eq.(^). The 
comparison yields an effective resistance R c s(uj) w 3g/ (4ir8(ui)). This indicates that the 
thermal currents actually flow in a skin layer right below to the metallic surface, in agreement 
with the estimates made by the NIST group [|lj. 

The ion heating rate obtained from the spectrum ([l2]) is plotted in figure 2. We observe that 

2 Rigorously speaking, the electric field also contains a contribution due to modes impinging from infinity in the 
vacuum space. In the short distance limit z <C A relevant here, however, this contribution may be shown to be very 
small (C. Henkel, K. Joulain, R. Carminati, J. -J. Greffet, in preparation). 
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Fig. Fig. 2. Heating rate for a trapped ion (mass 40 amu, charge q — e) trapped in a harmonic trap 
(frequency Q/2n = 1 MHz) above a copper surface at 300 K. Dots: exact evaluation of the electric near 
field fiuctations; solid line: asymptotic expansion (p~2|); dashed line: heating rate from the Johnson noise 
spectrum (M) with R(lo) = 1 f2. 



the electric near field fluctuates at a noise level much larger than the blackbody field (dotted 
line). It is also apparent that life times shorter than 1 s are to be expected when the trap size 
gets into the micrometre range. The figure also shows that the life times in current traps are 
not limited by the near field fluctuations discussed here, the observed heating rates being much 
larger. The NIST group proposed a model including fluctuating electric patch potentials that 
may explain this discrepancy, but little is known about the dielectric properties of these electric 
surface domains in the relevant frequency range [0]. In view of recent theoretical and experi- 
mental investigations [El ^, ^, [TTJ |l7| |l8|], a detailed understanding of ion heating processes is 
still lacking. 



2.3 Magnetic field noise spectrum 

To conclude this section, we turn to a different type of trap: a static inhomogeneous magnetic 
field that traps a paramagnetic atom at local minima of the magnetic field strength. Such traps 
are routinely used for evaporative cooling Jl9|], and miniature versions close to wires or surfaces 
have been proposed [^, |2(J and actually realised [ plj , p2| , ^3, 24, 2^]. Thermally fluctuating 



magnetic fields now may flip the atomic spin, leaving the atom possibly in an unbound potential 
(see figure 3). The spin flip rate is easily obtained from perturbation theory and relates to the 
spectral density of the magnetic field fluctuations: 

IW(r) =^$>|// a |/)</Mi)Sf (r;^) (13) 

Q, (3 
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Fig. Fig. 3. Fluctuating magnetic fields flip the spin of an atom in a magnetic trap and put it onto a 
non-trapping potential surface. 
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Fig. Fig. 4. Spin flip rate (loss rate) for a paramagnetic atom (magnetic moment /j, = /ib) in a magnetic 
trap above a copper surface at 300 K. The Larmor frequency takes the values ul/2tt = 1MHz (curve 
(a)) and 100 MHz (curve (b)). The rate obtained from the blackbody spectrum is much smaller (about 
10 13 s _1 for^ L /27r = 100 MHz). 

where ujl = /xS trap (r)/?i is the Larmor frequency. The magnetic field spectrum may be 
calculated using the theory outlined above, and one finds the following asymptotic result [^J: 

B [Z ' UJ) ~ 16^2 c 4 p (1 _ e-W^r) z + 353 (w) ) ^> 

For a spin 1/2 particle, we thus get the trap loss rate shown in figure 4. We note that in mi- 
crometre size traps, the life time is limited to less than a second due to near field fluctuations, 
which should be an observable effect in current experiments. 
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3 Decoherence in wave guides 

We now turn to the influence of fluctuating near fields on the quasi-free motion of atoms in 
linear or planar wave guides |Q |(| EH], pEj ]. The atomic matter wave is scattered from spatial 
inhomogeneities of the perturbing field, thus changing its momentum. We assume again a 
statistical description of the scattering process. The typical momentum transfer is thus of the 
order of %/i where I is the correlation length of the field. Energy is not conserved because the 
field is fluctuating, and the maximum energy transfer is roughly limited by fc^T. Note that this 
is typically much larger than the kinetic energies involved in cold atomic clouds. 

3.1 Transport equation 

We want to characterise the evolution of the single-particle spatial density matrix (or coherence 
function) 

/ o(r; S ) = <^(r + i S )V(r-|s)) (15) 

where the average (. . .) is taken over the spatial and temporal fluctuations of a perturbing po- 
tential V(r, t). It is useful to introduce the Wigner transform of the density matrix 

W(r,p) = J f± e ^ s / h p(r; S ) (16) 

that may be interpreted as a quasi-probability distribution in phase space. 

Using second-order perturbation theory, assuming gaussian statistics for the perturbing po- 
tential and doing a multiple-scale expansion of the Bethe-Salpeter equation for the coherence 
function, we get the following transport equation [^6|] : 

(d t + -p ■ V r + F cxt • V p ) W(r, p) = 

J d D p' S v (p' - p; E p , - E p ) [W(r, p') - W(r, p)] (17) 

where D — 1, 2 is the dimension of the wave guide and 5(q; AE) the spectral density of the 
fluctuating potential 

1 r r! D i (It 

S v (q; &E) = -^J J^d ( V ( r + t + T ) V ( r < *)> e- l(q - s - ABr)/R . (18) 

We have assumed that the potential is statistically stationary in both space and time. 

The left hand side of the transport equation ( |l7| ) describes the ballistic motion of the atom 
subject to the external (deterministic) force F cxt . The right hand side describes the scattering 
off the fluctuating potential. As a function of the momentum transfer q, e.g., the spectral density 
Sy (q; AE) is proportional to the spatial Fourier transform of the potential, as to be expected 
from the Born approximation for the scattering process p — > p' = p + q. The transport 
equation thus combines in a self-consistent way ballistic motion and scattering processes. 
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Fig. Fig. 5. Magnetic field spectrum vs. frequency at distance 1 /im above at copper surface at 300 K. 
The frequency u is expressed via the electromagnetic wave length A = 2-kc/uj. Note: a kinetic energy of 
10 /jK corresponds to a wave length 14 km. 
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Fig. Fig. 6. Spatial (normalised) correlation function of the thermal magnetic near field above a metallic 
surface at frequency lu/2tt — 30 MHz. The separation s gives the distance between two observations 
points at the same height z above the surface. The other parameters are identical to the previous figure 5. 
Dots: exact evaluation, solid lines: asymptotic expansions in the short-distance regime. For even lower 
frequencies, the correlation function is essentially unchanged. 

3.2 Example: magnetic perturbation 

For illustration purposes, we show in figure 5 the magnetic near field spectrum at a distance 
z = 1 /im above a flat metallic surface. This spectrum is proportional to Sy(q; AE) in ( |l7| ) 
for a planar magnetic waveguide. One observes that for typical kinetic energies of cold atoms, 
the magnetic spectrum is essentially flat. In the following, we shall hence approximate the 
perturbing field by a white noise. 

In figure 6, we show the normalised spatial correlation function of the thermal magnetic field 
above a metallic surface. One observes that in a planar wave guide above the surface, the field's 
spatial correlation length is of the order of the height z. The correlations decay algebraically 
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with the lateral separation s (measured in the waveguide plane parallel to the surface). 

3.3 Analytic solution of the transport equation 

For a broad band spectrum of the perturbation, we may neglect the dependence of SV(q, AE) 
on AE. The transport equation §17\ ) then simplifies because the integration over the scattered 
momentum p' is not restricted by energy conservation. Taking the Fourier transform with 
respect to both variables r and p (with conjugate variables k and s), it is simple to derive the 
following solution 



W(k, s; t) = Wo(k, s - Kkt/m) e - iF ««' s */ R 



x exp 



-7 f (1 - C(s - Kkt'/m))dt' 
Jo 



(19) 



Here, Wq (k, s) is the double Fourier transform of the Wigner function at initial time t = 0, and 
7 and the normalised spatial correlation function C (s) are related to the correlation function of 
the perturbation by 

7 C(s) = i / + dr (V(s, t) V(0, 0)), C(0) = 1. (20) 

n j_oo 

We also note that 7 is the scattering rate 7(p — * p') for 'forward scattering' processes where 
the final momentum p' approaches the initial p. For a magnetic wave guide above a metallic 
surface, the rate 7 is essentially of the same order of magnitude as the spin flip rate shown in 
figure 4. 

From the analytic solution (|l9|), it is easily checked that in the absence of the perturbation , 
the spatial width 8r 2 {t) of a cloud increases ballistically according to Sr 2 (t) = Sp 2 (0) t 2 /m 2 
where Sp(0) is the initial width of the cloud in momentum space (this latter width remains 
constant in this case, of course). 

3.4 Discussion 

Spatial decoherence. More interesting information may be obtained for a nonzero scattering 
rate 7. Note that the spatially averaged atomic coherence function is given by 

r(s) = Jd D rp(r;s) = WTk = 0,s) (21) 

The solution ( p"9[ ) therefore implies that the spatial coherence decays exponentially with time: 



r(s;i) = r (s)exp -yt(l - C(s)) - iF -st/h 



(22) 



The decoherence rate depends on the spatial separation between the points where the atomic 
wave function is probed, and is given by 7(s) = 7(1 — C(s)). It hence saturates to the value 7 
at large separations and decreases to zero for s — > 0. The decay of the coherence function ( |22| ) 
is illustrated in figure 7. One observes that at time scales t > 1/7, the spatial coherence is 
reduced to a coherence length £ coh ~ I. After a few collisions with the fluctuating magnetic 
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Fig. Fig. 7. Illustration of spatial decoherence in an atomic wave guide. The spatially averaged coherence 
function T(s, t) is plotted vs. the separation s for a few times t. Space is scaled to the field correlation 
length £ and time to the scattering time I/7. A Lorentzian correlation function for the perturbation is 
assumed. 



field, the long-scale coherence of the atomic wave function is thus lost and persists only over 
scales smaller than the field's correlation length (where different points of the wave function 
'see' essentially the same fluctuations). For larger times t I/7, decoherence proceeds at a 
smaller rate that is related to momentum diffusion, as we shall see now. 

Momentum diffusion at long times. The behaviour of the atomic momentum distribution at 
long times may be extracted from an expansion of the analytic solution (^2|) for small values of 
s. Assuming a quadratic dependence of the field's correlation function, C(s) w 1 — s 2 /£ 2 , as 
one would expect for Lorentzian correlations, we find that the atomic momentum distribution 
is gaussian at long times; it is centered at po + F cxt £ due to the external force, and its width 
increases according to a diffusion process in momentum space 

5p 2 (t)^S P 2 (0) + ^ (23) 

This was to be expected: the atoms perform a random walk in momentum space, exchanging a 
momentum of order Ti/£ per scattering time 1/7. The momentum diffusion coefficient D p = 
h 2 j/£ 2 that may be read off from ( |23| ) is consistent with this intuitive interpretation. Physically 
speaking, the atomic cloud is 'heated up' due to the scattering from the fluctuating potential. 
We note that the rate of change of the atomic kinetic energy in the wave guide plane is the 
same as the one for the tightly bound motion perpendicular to the metallic surface (see [^|] for 
a calculation of this rate). 

Translating the width of the momentum distribution into a spatial coherence length, we 
find a power-law decay at long times, £ co h = £/y/^ft. Finally, a similar calculation yields 
the width of the atomic cloud in position space: it increases 'super-ballistically' at long times, 
5r 2 (t) oc t 3 , as a consequence of heating. 
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4 Conclusion 

Particles in small traps close to macroscopic bodies are subject to fluctuating near fields that 
show noise spectra orders of magnitude above the blackbody level. This is because the ge- 
ometric distances involved are typically much smaller than the electromagnetic wave lengths 
associated with the relevant frequencies. As a consequence, the ground state of the vibrational 
motion of trapped ions is unstable, and coherences between different oscillator levels decay. We 
have developed a theoretical framework to compute the corresponding heating and decay rates. 
As a second consequence, quasi-free matter waves in a linear or planar wave guide in the vicin- 
ity of macroscopic bodies (substrates or wires) are scattered and lose their spatial coherence. 
We have identified the relevant time scale and obtained analytic estimates for the behaviour of 
the atomic coherence function at large spatial and temporal scales. These estimates should be 
useful, we hope, to design integrated atom optical circuits with controlled decoherence. 

Questions that could be addressed in the future pertain to detailed theories for trapped ion 
heating, as well as to transport processes for condensed atomic samples. Finally, the inclusion 
of interference effects in multiple scattering might provide a link to study weak and/or strong 
localisation of matter waves in wave guides. 

Acknowledgments. C.H. would like to thank R. Carminati and K. Joulain of the group of 
J. -J. Greffet at Ecole Centrale (Paris) for useful discussions. Furthermore, D. Leibfried, E. 
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